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1.0  BACKGROUND 


In  an  earlier  publication  (Bocker  &  Frieden,  1992)  a  new  covariant  matrix 
representation  of  classical  electromagnetic  theory  for  vacuum  was  presented.  The 
quintessential  basis  of  this  representation  is  a  skew-Hermitian  space-time  8-by-8 
differential  matrix  operator 
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and 


x  =  ict, 


(1) 


(2) 


(3) 


The  symbol  c  represents  the  speed  of  light  in  vacuum.  Use  of  the  square  matrix  operator 
[M]  allows  Maxwell's  equations  to  be  placed  in  a  compact  matrix  form.  From  the  matrix 
form  of  Maxwell's  equations,  other  landmark  effects  of  electromagnetic  theory  are  easily 
derived  with  use  of  the  simple  matrix  multiply  operation:  (a)  the  electromagnetic  wave 
and  charge  continuity  equations;  (b)  the  Lorentz-gauge  and  Coulomb-gauge  definitions 
of  the  electromagnetic  potentials;  (c)  the  wave  equations  for  the  potentials;  and  (d) 
Poynting's  theorem  on  energy  conservation.  Taking  the  four-dimensional  Fourier 
transform  of  the  matrix  form  of  Maxwell's  equations  leads  to:  (e)  a  Fourier 
representation  of  Maxwell's  equations;  (f)  their  inversion,  for  the  fields  directly  in 
terms  of  the  sources  in  Fourier-space;  and  (g)  corresponding  inversion  formulae  in 
direct-space  through  the  use  of  the  convolution  theorem. 
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Because  of  the  power  of  matrix  operations,  we  were  able  to  derive  the  key  effects  of 
electromagnetic  theory  without  the  need  for  the  usual  plethora  of  vector  calculus 
identities  that  have  become  the  standard  in  these  derivations:  Stokes'  theorem,  the 
Divergence  theorem.  Green's  theorem,  and  the  formula  for  the  curl  of  a  curl,  etc.  Instead, 
the  simple  matrix  operations  of  matrix  multiplication  and  matrix  inversion  were  used. 

In  this  document,  the  matrix  representation  of  classical  electromagnetic  theory  for 
vacuum  will  be  extended  to  include  the  presence  of  matter.  Emphasis  is  placed  on 
electromagnetic  wave  propagation  in  linear,  homogeneous,  anisotropic  optical  media  of 
infinite  extent  v'^hout  boundaries.  A  subsequent  document  treating  wave  propagation 
in  optical  media  will  include  boundaries. 

2.0  EXTENSION  OF  THE  MATRIX  FORMULATION 
2.1  MAXWELL  FIELD  EQUATIONS 

The  fundamental  equations  of  classical  electromagnetic  phenomena,  namely  the 
Maxwell  field  equations,  serve  as  our  starting  point.  In  the  Gaussian  system  of  units,  the 
four  Maxwell  field  equations  in  vector  form  are  given  by  the  following  (Jackson,  1962): 

Ampere-Maxwell  law 

VxH(r.f)  =  i|-D(r,f)+—  f(r,t),  (4a) 

Cdt  c 

Gauss'  law  for  electricity 

V*D(r,f)  =4npe(r,t),  (4b) 

Faraday's  law  of  induction 

Vx  E  (r,  t)  =  -  -|-B  (r,  t)  -  —Jm  ( r ,  t ) ,  (4c) 

Cdt  C 

and 

Gauss'  law  for  magnetism 

V«B  (r,  t)  =  4rrp m(r,t).  (4d) 
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In  rectangular  coordinates. 


V 


(A  A  A) 

W  dy  dz} 


(5) 


and 


r  =  (x,  y,  z) .  (6) 

The  physical  quantities  appearing  in  the  Maxwell  field  equations  are  E{r,t )  the 
electric  field;  D  ( r ,  t )  the  electric  displacement;  B  ( r ,  t )  the  magnetic  induction;  H  ( r,t ) 
the  magnetic  field;  f  ( r ,  f)  the  electric  current  density;  Jm  (r,  t )  the  magnetic  current 
density;  p e  (r,  t )  the  electric  charge  density;  pm  (r,  t)  the  magnetic  charge  density; 
and  (r,t)  a  space- time  point.  Both  magnetic  charge  and  current  densities  (Magid,  1972) 
have  been  included  in  Maxwell's  equations  for  purposes  of  completeness.  They,  of 
course,  can  be  set  equal  to  zero  since  magnetic  charge  has  not  been  discovered  in  nature. 

The  electric  displacement  and  electric  field,  as  well  as  the  magnetic  induction  and 
magnetic  field,  are  related  (Jackson,  1962)  through  the  expressions 


D(r,t)  =  E(r,t)+4n  P(r,t)  (7a) 

and 

B(r,t )  =  H(r,t) +4n  M(r,t) ,  (7b) 

where  P(r,t )  is  the  macroscopic  polarization  and  M(r,t)  is  the  macroscopic 
magnetization.  By  using  the  approach  adopted  previously  by  Bocker  and  Frieden  (1992), 
the  four  Maxwell  field  equations  (4)  with  the  use  of  equation  (7)  can  be  cast  into  the 
following  matrix  form 
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The  4-by-4  matrix  operators 


M. 


and 


represents  the  4-by-4  null  matrix.  The 


M2j  are  defined  by  equation  (2).  The  matrix  [O] 
4-by-l  field  vectors  appearing  in  equation  (8)  are 
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defined  by 
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The  source  vectors  in  equation  (8)  are  ioativistic  4-vectors  defined  by 
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2.2  ELECTROMAGNETIC  FIELD  WAVE  AND  CHARGE 
CONTINUITY  EQUATIONS 

Multiply  both  sides  of  matrix  equation  (8)  by  the  complex  conjugate  of  the  space-time 
operator  [M].  This  gives 
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+  471 

Di  D2 

O 

D 

A 

d2  d1 

A 

4k 
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The  matrix  [D]  is  defined  by 


[D]  ee  +  [M2\*[M2 J 


□  2  0  0  0 

0  D2  0  0 

0  0  D2  0 

0  0  0  D2 


(11) 


(12) 
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The  D'Alembertian  operator  D2  appearing  in  equation  (12)  is  defined  by  (Ohanian, 
1988) 


□2 


B  a2  a2  a2_  a2 

dx 2  dy2  dz2  3t2 


The  matrices  [D^]  and  [D2]  are  defined  by 


(13) 


[D,]  s  [MjAm,]  and  [D2]  s  [M2]  *  [M^  . 


(14) 


Again,  the  matrix  [O]  represents  the  4-by-4  null  matrix.  Matrix  equation  (11)  is 
equivalent  to  eight  scalar  equations.  These  scalar  equations  are  equivalent  to  the  vector 
electromagnetic  field  wave  equations: 


VZE  - 


C2dt2 


a2E  =  4rc  Vpe+~!^+^Vx/m 


4k  ^ 2  4k 

-4k  V  (V*P)  +  — —rP  + 


c2dt2'  c3t 


Ir(VxM) 


(15a) 
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1  32  ii- 


c2dt2 


H  =  4k  Vpm  +  4J|r/w-^Vx/f 


,2  at7 


4k  4K  3 

-4k  V  (V*M)  +  — ^-(VxP) , 

c 2  dt‘ 


4k 


.2-w2‘"  c  dt 


and  the  electric  and  magnetic  charge  continuity  equations 
+  pc  =  0 

and 

V«/w  +  |-pm  =  0. 

J  dV 


(15b) 


(16a) 


(16b) 
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2.3  ELECTROMAGNETIC  POTENTIALS  AND  LORENTZ 
CONDITIONS 

We  found  before  that  the  space-time  operator  [M]  defines  the  eight  scalar  Maxwell 
field  equations.  We  now  observe  the  complex  conjugate  of  [M]  also  provides  the 
definition  of  the  electromagnetic  fields  in  terms  of  the  familiar  vector  and  scalar 
potentials.  In  particular,  for  the  Lorentz  gauge  we  have 

A 
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where  the  relativistic  4-vector  potentials  are  defined  by 
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(18) 


Equation  (17)  is  equivalent  to  eight  scalar  equations.  Six  of  these  are  equivalent  to  the 
two  vector  equations  defining  the  relationship  between  the  electromagnetic  potentials 
and  the  electromagnetic  fields: 


6  1  d  A  &  T7.  .  A 


and 


E  =  -  V  <p  — ~-A  -  Vx  A 
c  at 

H  =  -  V  <pm  -  -|-Am  +  Vx  Ae. 
C  at 


The  remaining  two  scalar  equations  correspond  to  the  Lorentz  conditions 


and 
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V'Am  +  -^-<pm  =  0. 
C  at 


(19a) 

(19b) 

(20a) 

(20b) 


Equations  (19)  and  (20)  are  easily  verified  by  explicit  multiplication  of  the  right-hand 
side  of  matrix  equation  (17). 


2.4  ELECTROMAGNETIC  POTENTIAL  WAVE  EQUATIONS 


We  know  the  electromagnetic  vector  and  scalar  potentials  satisfy  inhomogeneous 
wave  equations  (Jackson,  1962).  Once  again,  we  show  this  by  simple  matrix 
multiplication.  First,  substitute  the  matrix  expression  (equation  17)  into  the  matrix 
representation  (equation  8)  of  the  Maxwell  field  equations.  This  gives 
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For  the  Lorentz  gauge  choice  of  [M], 
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by  direct  evaluation  of  the  matrix  product.  Then,  by  equations  (21)  and  (22),  the 
following  matrix  representation  of  the  electromagnetic  potential  wave  equations  is 
obtained 
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By  explicit  evaluation,  this  single  matrix  equation  is  equivalent  to  eight  scalar 
equations.  Six  of  these  scalar  equations  are  equivalent  to  the  following  pair  of 
inhomogeneous  wave  equations  for  the  electric  and  magnetic  vector  potentials: 
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The  remaining  two  represent  the  inhomogeneous  wave  equations  for  the  electric  and 
magnetic  scalar  potentials: 


V2cpe  -  ^r-^r-(pe  =  -4  n  pe  +  4n  V»P 
cdt 2 

and 

V2<pm-I-3i<pm  =  -4*  pm  +  4n  V.M. 
c2S  l2 


(25a) 


(25b) 


2.5  ENERGY  CONSERVATION:  POYNTING'S  THEOREM 

The  law  of  conservation  of  energy,  often  called  Poynting's  theorem  (Jackson,  1962),  is 
an  important  milestone  of  electromagnetic  theory.  Again,  simple  matrix  manipulation  of 
equation  (8)  will  accomplish  the  derivation.  Multiply  both  sides  of  equation  (8)  by  the 
Hermitian  conjugate  of  the  electromagnetic  field  vector.  Directly, 


A 
A 

This  is  a  single  scalar  equation  representing  Poynting's  theorem 


Ml  M2 
M2  Ml 


+  4n 


M1  O 
O  M, 


4k 
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V*  S  +  Jp-w  +  E  •  f  +  H  •  Jm  +  E  •^  +  H  •  =  0. 

dt  ’  ’  dt  dt 


dM 


(26) 


(27) 


We  define  quantity  u ,  the  total  energy  density  of  the  electromagnetic  fields,  to  obey 

u  =  ~(E»E  +  H»H)  (28) 

8k 

and  define  S ,  the  Poynting  vector  representing  energy  flow,  by 


S=^(E  x  H) .  (29) 

471 

Again,  this  is  verified  by  explicit  evaluation  of  the  matrix  products  in  equation  (26). 
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3.0  OPTICAL  MATERIALS 


The  study  of  the  propagation  of  electromagnetic  radiation  through  matter  comprises 
an  important  branch  of  modern  optics  (Fowles,  1968).  Many  of  the  optical  properties  of 
matter  can  be  understood  on  the  basis  of  classical  electrodynamics.  In  this  section,  the 
matrix  representations  of  both  the  Maxwell  field  equations  (8)  and  electromagnetic  wave 
and  charge  continuity  equations  (11)  will  be  used  to  mathematically  describe  wave 
propagation  in  linear,  homogeneous,  anisotropic  optical  media  of  infinite  extent. 
Monochromatic  plane-wave  solutions  in  the  absence  of  charges  and  currents  will  be 
considered. 

3.1  MAXWELL  FIELD  EQUATIONS 

For  linear,  homogeneous,  anisotropic  optical  media,  both  the  macroscopic 
polarization  and  the  magnetization  vectors  are  related  to  the  electric  and  the  magnetic 
field  vectors  through  4-by-4  electric  and  magnetic  susceptibility  tensors  [x£]  and  [x  ]  • 
Mathematically,  we  can  express  these  relationships  in  the  compact  matrix  form 
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In  a  medium  with  neither  free  charges  nor  currents,  we  also  have 


(30) 


(31) 


Substitution  of  equations  (30)  and  (31)  into  the  Maxwell  field  equations  (8)  gives 
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For  monochromatic  plane-wave  fields,  the  electromagnetic  field  vector  in  equation 
(32)  can  be  expressed  in  the  form 
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where  the  relativistic  4-vectors  [r]  and  [k]  are  defined  by 
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(34) 


The  superscript,  T,  appearing  in  equation  (33)  denotes  transpose.  The  4-by-l  vectors 
[/ol]  and  [/o2]  appearing  in  equation  (33)  are  both  constant  and  describe  the 
polarization  properties  of  the  plane-wave  field.  With  the  use  of  equation  (34),  the 
argument  of  the  exponential  function  in  equation  (33)  is  given  by 


[k]T[r] 


(35) 


The  wave  vector  k  in  equation  (35)  defines  the  direction  of  wave  propagation  within  the 
medium.  Its  magnitude  k,  commonly  referred  to  as  the  wavenumber,  is  related  to  the 
wavelength  X  of  the  electromagnetic  wave  through  the  equation 

k  =  (36) 

The  angular  frequency  o)  appearing  in  equation  (34)  is  related  to  the  frequency  / 
through  th°  equation 

co  =  2  ft/.  (37) 


Note  the  monochromatic  plane-wave  representation  (equation  33)  is  equivalent  to  the 
vector  equations 
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and  H(r,t )  = 
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(38) 


Substitution  of  equation  (33)  into  equation  (32)  gives  the  following  monochromatic 
plane-wave  matrix  representation  of  the  Maxwell  field  equations  for  an  infinite 
anisotropic  optical  medium: 
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The  4-by-4  matrices 


K, 


and 


K. 


are  defined  by 
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where,  by  definition, 


(41) 


Note  the  matrix  equation  (39)  represents  a  system  of  algebraic  equations  in  contrast  to 
equation  (32),  which  describes  a  system  of  first  order  partial  differential  equations. 
Through  the  use  of  matrix  multiplication,  equation  (39)  can  be  rewritten  in  the  form 
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(42) 


The  4-by-4  dielectric  tensor  [e]  is  related  to  the  electric  susceptibility  tensor  [ xg] 
and  the  4-by-4  identity  matrix  [I]  through  the  defining  equation 


[e]  s  [I]  +4n[xey 


(43) 


Similarly,  the  4-by-4  permeability  tensor  [p]  is  related  to  the  magnetic  susceptibility 
tensor  fxml  and  the  identity  matrix  [I]  through  the  defining  equation 

[p]  s  [I]  +47i^%mJ.  (44) 

Using  definitions  (43)  and  (44)  with  equation  (42)  gives 
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Matrix  equation  (45)  is  equivalent  to  the  following  transversality  conditions  in  vector 
form: 


and 


kxH  = — — D  , 

(46a) 
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o 

(46b) 
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k  x  £  =  +  -  B  , 
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(46c) 
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0 

(46d) 

The  next  important  step  in  this  development  is  to  cast  equation  (45)  into  an 
eigenvalue  equation.  The  4-by~4  matrices  [K^]  and  [K2] ,  first  defined  in  equation  (40), 
can  be  rewritten  in  the  form 

co 

[KJ  =  iTclaj]  +-  [I]  and  [K2]  =  ik  [a2] ,  (47) 


where  the  matrices  [o^]  and  [a2]  are  defined  by 
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Again,  the  matrix  [I]  represents  the  4-by-4  identity  matrix.  The  matrix  elements 
a  ,  a  ,  and  a  appearing  in  equation  (48)  are  direction  cosines  (Morrill,  1961)  defining 

**  p  #  ^ 

the  direction  of  the  wavevector  k.  These  direction  cosines  are  related  to  the  components 
of  k  and  the  wavenumber  k  through  the  equation 


k  =  ka  , 
x  x 

Also, 


and  k  =  ka  . 
2  2 


k2  =  k2  +  k2+k2  and 
x  y  z 


111 
a  +a'  +a  =  1. 
x  y  z 


(49) 


(50) 
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With  the  use  of  equations  (47)  and  (48),  matrix  equation  (45)  can  be  rewritten  as  a 
generalized  eigenvalue  equation 


c^e 

a2 

4, 

i 

e 

0 

4," 

a2 

04  B 

[o2 

n 

0 

ft. 

^02 

where  the  index  of  refraction  n  is  related  to  the  temporal  frequency  /,  the  wavelength  X, 
and  the  speed  of  light  in  vacuum  c  through  the  equation 

nXf  =  c.  (52) 

The  matrix  eigenvalue  representation  (equation  51)  of  the  Maxwell  field  equations, 
hereafter  referred  to  as  the  Maxwell  eigenvalue  matrix  representation,  is  the  most 
important  result  of  subsection  3.1.  Knowledge  of  both  the  dielectric  and  permeability 
tensors  [e]  and  [p]  of  the  optical  medium  of  interest,  as  well  as  the  direction  of  wave 
propagation  within  the  optical  medium,  completely  specifies  the  two  8-by-8  matrices 
appearing  on  both  sides  of  equation  (51).  Once  these  8-by-8  matrices  have  been  specified, 
then  the  unknown  eigenvalues  i/n  (hence,  eigenindices  of  refraction  n)  and 
corresponding  unknown  eigenvectors  [/  /  ]  7  (hence,  eigenstates  of  polarization) 

can,  in  principle,  be  determined.  Off-the-shelf  computer  software  packages,  like 
MATLAB  for  numerical  computations  and  MATHEMATICA  for  symbolic 
manipulations,  can  facilitate  the  process  of  finding  the  eigenvalues  and  eigenvectors  of 
equation  (51). 


3.2  ELECTROMAGNETIC  FIELD  WAVE  AND  CHARGE 
CONTINUITY  EQUATIONS 


As  previously  demonstrated,  multiplication  of  both  sides  of  the  matrix  form  of  the 
Maxwell  field  equations  (8)  by  the  complex  conjugate  of  the  space- time  operator  [M] 
leads  to  the  matrix  form  of  the  electromagnetic  field  wave  and  charge  continuity 
equation  (11).  For  the  case  of  monochromatic  plane-wave  fields  (equation  33) 
propagating  in  linear,  homogeneous,  anisotropic  optical  media  (equation  30)  in  the 
absence  of  charges  and  currents  (equation  31),  the  matrix  equation  (11)  simplifies  to 


* 

r 

k2 

KjE  K2 

4i 

0 

K2 

h 

k2 

l°2. 

0 

(53) 
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Matrix  multiplication  further  yields 


*  * 

*  * 

r  -i 

r  -i 

iC-j  X|£  +  K  2  K  2 
*  * 

K^K2  +  K2K^  p 

*  * 

4i 

[o  2 

= 

0 

K2K1e  +  K1K2 

k2k2  +  k^k^ 

0 

Substituting  the  mathematical  expressions  for  [K-j]  and  [/(.,]  from  equation  (47)  into 
equation  (54),  with  the  help  of  equation  (52),  leads  to  the  following  important  result  for 
subsection  3.2: 


I-a  a  (e-I)  -a  (p-I) 

11  n  2 

-a  (e-I)  I-a1a1(p-I) 

n  2  11 

Hereafter,  equation  (55)  will  be  referred  to  as  the  electromagnetic  wave  equation 
matrix  representation.  Equation  (55)  greatly  simplifies  in  structure  when  either  the 
dielectric  tensor  [e]  or  the  permeability  tensor  [|i]  is  equal  to  the  identity  matrix  [I], 
This  fact  will  become  quite  evident  in  section  4  where  monochromatic  plane-wave 
propagation  in  a  variety  of  nonmagnetic  optical  media  is  considered  in  greater  detail. 

4.0  SPECIFIC  APPLICATIONS 

The  utility  of  using  the  Maxwell  eigenvalue  matrix  representation  (equation  51)  or 
the  electromagnetic  wave  equation  matrix  representation  (equation  55)  in  solving  wave- 
propagation  problems  for  a  variety  of  nonmagnetic  optical  media  will  be  pursued  in  this 
section.  Examples  of  nonmagnetic  optical  materials  considered  include  crystalline 
materials,  optically  active  materials,  and  electrooptic  materials.  For  nonmagnetic  optical 
materials,  the  permeability  tensor,  defined  by  equation  (44),  simplifies  to  the  following 
form 


Ol 


02 


n 


e  O 
O  p 


oi 


02 


(55) 


[p]  =  [I].  (56) 

At  this  point,  we  introduce  the  optical  impermeability  tensor  [  rj  ]  (Yariv  &  Yeh,  1984), 
which  is  defined  as  the  multiplicative  inverse  of  the  dielectric  tensor  [e]  .  That  is, 

[TtfMer1.  (57) 
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With  the  use  of  equations  (56)  and  (57),  the  Maxwell  eigenvalue  matrix  representation 
(equation  51)  can  be  rewritten  in  the  form 


ric^e 

ia2 

u 

i 

4i 

a2 

v 

fo2 

n 

f()2 

(58) 


Matrix  multiplication  involving  the  left-hand  side  of  equation  (58)  leads  to  the  following 
pair  of  matrix  equations: 

+  !«,)  I/o2)  -  X1  (59a) 

and 

til  [«,1  [el  [/„,]  +  [ril  [a,]  [/o.]  =  i  [/o]l .  (59b) 


Similarly,  with  the  use  of  equations  (56)  and  (57),  the  electromagnetic  wave  equation 
matrix  representation  (equation  55)  simplifies  to 

y\ -Tia1a1  (e- 1)  O 

(i/n)  a2(e-I)  I 

Matrix  multiplication  involving  the  left-hand  side  of  equation  60  leads  to  the  following 
pair  of  matrix  equations: 


4i 

1 

4, 

fo  2 

=  72 

fo2 

(60) 


M  <[I] -[<*,][<*, ne-Il>[/0])  =^[/0]J  (61a) 

and 

[f02]  =-T^7[a2lle-I][/01].  (615) 

n  - 1 

Equation  (61a),  like  equation  (58),  is  an  eigenvalue  equation.  However,  equation  (61a) 
involves  the  unknown  vector  [/ol  j  only.  Equation  (61b)  allows  [/o-J  to  be  determined 
directly  in  terms  of  the  solution  f/ol]  obtained  from  equation  (61a).  Thus,  for  non¬ 
magnetic  materials,  the  electromagnetic  wave  equation  matrix  representation  (equation 
55)  reduces  to  the  simplify  pair  of  matrix  equations  given  by  equation  (61). 
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At  this  time,  it  is  important  to  emphasize  the  following  point.  Either  the  Maxwell 
eigenvalue  matrix  representation  (equation  58  or  59)  or  the  electromagnetic  wave 
equation  matrix  representation  (equation  60  or  61)  can  be  used  to  determine  the 
eigenvalues  (hence,  eigenindices  of  refraction)  and  corresponding  eigenvectors  (hence, 
eigenstates  of  polarization)  associated  with  the  propagation  of  monochromatic 
electromagnetic  plane-wave  radiation  in  nonmagnetic  media. 


4.1  CRYSTALLINE  MATERIALS 

A  basic  feature  of  a  crysf  lline  material,  as  far  as  optical  properties  are  concerned,  is 
that  crystals  are  generally  electrically  anisotropic.  Hence,  the  macroscopic  polarization 
produced  in  the  crystal  by  an  applied  electric  field  varies  in  a  manner  that  depends  on 
the  direction  of  the  applied  electric  field  in  relation  to  the  crystal  lattice.  It  is  a  simple 
matter  to  show  that  for  ordinary  nonabsorbing  crystals  (Born  &  Wolf,  1965),  there 
always  exists  a  set  of  coordinate  axes,  called  principal  axes,  such  that  the  dielectric  tensor 
[e]  assumes  the  diagonal  form 


[e] 


£n  0  0  0 

0  e22  0  0 

0  0  e330 

0  0  0  1 


(62) 


Using  equation  (57),  we  find  the  optical  impermeability  tensor  is  given  by 


mi 


e'J  000 
0  0  0 

0  0  e33° 

0  0  0  1 


(63) 
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Substituting  these  expressions  for  [e]  and  [q]  back  into  the  Maxwell  eigenvalue  matrix 
representation  (equation  58)  yields  the  following  8-by-8  eigenvalue  equation: 


0 

0 

0 

-a  e1? 

X  11 

0 

-a  eJ 
2  ii 

a  e.3 

y  11 

0 

iE 

X 

iE 

X 

0 

0 

0 

-a  cl 
y  22 

a  e~l 

2  22 

0 

-a  e~] 
a:  22 

0 

iE 

y 

iE 

y 

0 

0 

0 

-a  el] 
2  33 

-a  el] 

y  33 

a  e"] 
a:  33 

0 

0 

iE 

2 

iE 

2 

a  e^ 
x  n 

y  22 

a  e,. 
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0 

0 

0 

0 

0 

i 

0 

0 

-a 

2 

a 

y 

0 

0 

0 

0 

-a 

X 

H 

X 

n 

H 

X 

a 

2 

0 

-a 

X 

0 

0 

0 

0 

-a 

y 

H 

y 

H 

y 

-a 

y 

a 

X 

0 

0 

0 

0 

0 

-a 

2 

H 

2 

H 

2 

0 

0 

0 

0 

a 

X 

a 

y 

a 

2 

0 

0 

0 

(64) 


Substituting  the  expressions  for  [e]  and  [T]]  back  into  the  electromagnetic  wave 
equation  matrix  representation  (equations  61a  and  61b)  yields  the  following  3-by-3 
eigenvalue  equation: 


l+a2(e11-l)  a  a  (e0_  —  1 )  a  a  (e__-l) 

X  11  X  y  22  2  X  33 

r 

r 

r 

ell  £11  ell 

L 

X 

E 

X 

a  a  (e^-1)  l  +  a2(e„.-l)  a  a  (e_,-l) 
x  y  ll  y  22  y  z  33 

E 

1 

E 

£22  £22  £22 

y 

/ 

n 

I—* 

y 

°zax(ell"1)  ay  <£22  -  ^ '  1  +  “z<e33-1) 

E 

*T 

E 

•J 

e33  e33  e33  J 

(65a) 
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as  well  as  the  solution  for  the  magnetic  field  in  terms  of  the  electric  field 


0  ~az  (e22  -  1) 

a,  (£jj  - 1)  0 

-ay(en-l)  a;c(e22-l) 


a  (e„ 
y  33 


0 


(65b) 


A  crystalline  medium  in  which  all  principal  dielectric  tensor  elements  are  different 
has  two  optic  axes  (Fowles,  V.  3).  For  this  case,  the  crystal  is  said  to  be  biaxial.  On  the 
other  hand,  two  of  the  principal  dielectric  tensor  elements  of  uniaxial  crystal  are  equal.  A 
uniaxial  crystal  has  a  single  optic  axis.  A  linear,  homogeneous,  isotropic  optical  medium 
is  commonly  referred  to  as  a  Class  A  dielectric.  The  principal  dielectric  tensor  elements 
of  a  dielectric  are  all  equal.  Vacuum  is  a  special  case  of  a  dielectric.  For  vacuum,  the 
principal  dielectric  tensor  elements  are  unity.  To  illustrate  the  use  of  equations  (64)  and 
(65),  let  us  consider  wave  propagation  in  a  vacuum,  a  dielectric,  and  a  uniaxial 
crystalline  medium,  respectively. 

4.1.1  Vacuum 

For  a  vacuous  medium,  the  elements  of  the  dielectric  tensor  [e]  are  given  by 

611  =  e22  =  e33  =  l'  (66) 

Substituting  these  values  of  en,  e22,  and  e33  back  into  equation  (65a)  gives 


The  3-by-3  matrix  on  the  left-side  of  equation  (67a)  is  independent  of  the  direction 
cosines.  Inspection  of  equation  (67)  indicates  the  only  physically  acceptable  value  for  the 
index  of  refraction  is 

n  =  1  (68) 

for  all  possible  states  of  polarization. 
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Substituting  the  values  of  en,  e22,  e33  and  n  from  equations  (66)  and  (68), 
respectively,  into  the  Maxwell  eigenvalue  matrix  representation  (equation  64)  gives 


0  0  0  -a  0  -a 

x  z 

0  0  0  -a  a  0 

y  2 

0  0  0  -a  -a  a 

z  y  x 

a  a  a  0  0  0 

x  y  z 

0  -a  a  0  0  0 

2  y 

a  0  -a  0  0  0 

2  X 

-a  a  0  0  0  0 

y  * 

0  0  0  0  a  a 

x  y 


(69) 


Matrix  equation  (69)  is  equivalent  to  eight  scalar  equations.  Six  of  these  scalar  equations 
correspond  to  the  vector  equations 


Eo  =  -(axHo)  (70a) 

and 

Ho  =  +(axEo),  (70b) 

and  the  remaining  two  scalar  equations  are  equivalent  to 

a  •  E  =  0 
o 

and 

a*H  =0. 
o 

a  is  the  unit  vector  whose  components  are  a^.,  a  ,  and  a7  .  Equations  (70)  and  (71) 
reveal  the  three  vectors  E°,  Hq,  and  a  are  mutually  orthogonal  to  one  another.  Also, 
the  magnitudes  of  the  electric  and  magnetic  field  vectors  are  equal.  That  is 


(71a) 

(71b) 


Thus,  monochromatic  electromagnetic  plane-waves  in  vacuum  propagate  at  the 
speed  of  light,  namely  c.  The  electric  and  magnetic  field  vectors  associated  with  these 
waves  are  orthogonal  to  each  other  as  well  as  to  the  direction  of  wave  propagation.  Such 
waves  are  called  transverse  waves.  Any  state  of  polarization  is  allowed  as  long  as  it 
satisfies  the  above  transversality  conditions. 
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4.1.2  Dielectrics 


As  previously  mentioned,  the  principal  dielectric  tensor  elements  of  a  dielectric 
material  are  all  equal 


e 


11 


£  =  £ 
t22  33 


£>  1. 


(73) 


Substituting  these  values  of  the  principal  dielectric  tensor  elements  back  into  the  3-by-3 
eigenvalue  equation  (65a)  gives 
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(74) 


The  only  physically  acceptable  value  for  the  index  of  refraction  satisfying  equation  (75)  is 

n  =  Je.  (75) 

This  value  of  the  index  of  refraction  applies  to  all  possible  states  of  polarization. 
Substituting  the  values  of  e^,  e22,  and  e33  back  into  the  8-by-8  eigenvalue  equation 
(64)  gives 
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Matrix  equation  (76)  is  equivalent  to  eight  scalar  equations.  Six  of  these  scalar  equations 
correspond  to  the  vector  equations 


and 


E  =4(axH  ) 

°  Ji  ° 


H  =  +Je(axE  ) , 
o  o 


and  the  remaining  two  scalar  equations  are  equivalent  to 


(77a) 

(77b) 


a*E  =  0 
o 

and 


(78a) 


a*H  =  0. 

o 


(78b) 


Equations  (77)  and  (78)  tell  us  the  three  vectors  E°,  H^,  and  a  are  mutually  orthogonal 
to  one  another.  Also,  the  magnitudes  of  the  electric  and  magnetic  field  vectors  are  related 
by  the  equation 


H  =  Je  E  .  (79) 

o  v  o 

So  monochromatic  electromagnetic  plane-waves  in  a  dielectric  optical  medium 
propagate  at  the  speed  of  c/n,  where  the  index  of  refraction  n  is  given  by  equation  (75). 
The  electric  and  magnetic  field  vectors  associated  with  these  waves  are  orthogonal  to 
each  other  as  well  as  the  direction  of  wave  propagation.  Any  state  of  polarization  is 
allowed  providing  it  satisfies  the  transversality  conditions  (equations  77  and  78). 


4.1.3  Uniaxial  Crystals 


As  previously  indicated,  two  of  the  principal  dielectric  tensor  elements  of  a  uniaxial 
crystal  are  equal.  Also,  a  uniaxial  crystal  has  a  single  optic  axis.  Without  loss  of 
generality,  we  will  take  the  z-axis  as  the  optic  axis  of  the  uniaxial  crystal.  The  dielectric 
tensor  elements  then  have  the  form 


£11  E22  £o 


and 


e,,  =  e  . 
33  e 


(80) 


is  referred  to  as  the  ordinary  dielectric  constant  and  e  is  called  the  extraordinary 
dielectric  constant. 
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For  both  uniaxial  and  biaxial  crystals,  the  speed  of  propagation  of  a  light  wave  in  the 
crystal  is  a  function  of  the  direction  of  propagation  and  the  polarization  of  the  light. 
There  exists,  in  general,  two  possible  values  of  the  phase  velocity  for  a  given  direction  of 
propagation.  These  two  values  are  associated  with  mutually  orthogonal  polarizations  of 
the  light  waves.  Crystals  are  thus  said  to  be  doubly  refracting  or  birefringent.  The  8-by-8 
eigenvalue  equation  (64)  or  the  3-by-3  eigenvalue  equation  (65a),  in  conjunction  with 
equation  (65b),  can  be  used  to  find  these  phase  velocities  and  corresponding 
polarizations  states.  Let  us  now  look  at  some  specific  examples  to  ilb'strate  these  ideas. 


Example  1:  Wave  Propagation  Along  the  x-Axis.  The  direction  cosines  describing  wave 
propagation  in  the  positive  x-direction  are  given  by 


ax  =  h  ay  =  °'  “z  =  °‘  (81) 

Substituting  the  dielectric  tensor  elements  from  equation  (80)  and  the  direction  cosines 
from  equation  (81)  back  into  the  8-by-8  eigenvalue  equation  (64)  gives 
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(82) 


It  is  a  simple  matter  to  show  that  only  two  physically  acceptable  eigenvalues  satisfy 
equation  (82).  The  corresponding  eigenindices  of  refraction  are  given  by 
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Associated  with  each  value  of  n  is  a  corresponding  eigenvector  solution  (eigenstate  of 
polarization).  The  eigenvector  solutions  are 
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where  C  is  an  arbitrary  cons’  .nt  These  two  eigenvector  solution  present  orthogonal 
states  of  linearly  polarized  'iectron.agnetic  radiation 

Substituting  the  dielectric  tensor  elements  from  equation  (80)  and  the  direction 
cosines  from  equation  (81)  bad  nto  equations  (65a)  and  (65b),  respectively,  gives 
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It  can  be  shown  that  only  two  physically  acceptable  eigenvalues  satisfy  equation  (85a). 
The  corresponding  eigenindices  of  refraction  are  again  given  by  equation  (83).  The 
corresponding  eigenvectors  (electric  fields)  obtained  from  equation  (85a),  and 
subsequent  solutions  (magnetic  fields)  obtained  from  equation  (85b),  are  identical  to  the 
electric  and  magnetic  field  solutions  summarized  in  equation  (84).  We  see  the  Maxwell 
eigenvalue  matrix  representation  (equation  64)  and  the  electromagnetic  wave  equation 
matrix  representation  (equation  65)  yield  identical  results.  As  previously  indicated, 
either  representation  can  be  used  in  solving  wave-propagation  problems  in  crystalline 
media. 


Example  2:  Wave  Propagation  in  the  xz-Plane.  Our  second  example  concerns  wave 
propagation  in  the  xz-plane.  Suppose  the  direction  of  wave  propagation  is  such  that  the 
wavevector  k  makes  an  angle  of  45  degrees  with  respect  to  both  the  positive  x-  and  z- 
axes.  The  direction  cosines  associated  with  this  example  are  given  by 

Jl  „  J2 

%=T'  V0'  “z=T  (8« 

Substituting  the  dielectric  tensor  elements  from  equation  (80)  and  the  direction  cosines 
from  equation  (86)  back  into  the  8-by-8  eigenvalue  equation  (64)  gives 
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It  can  be  easily  shown  that  only  two  physically  acceptable  eigenvalues  satisfy  equation 
(87).  The  corresponding  eigenindices  of  refraction  are  given  by  the  following  formula: 
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Associated  with  each  value  of  n  is  a  corresponding  eigenvector  solution  (eigenstate  of 
polarization).  The  eigenvector  solutions  are 
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where  C  is  an  arbitrary  constant.  These  two  eigenvector  solutions  represent  orthogonal 
states  of  linearly  polarized  electromagnetic  radiation. 


Example  3:  Wave  Propagation  Along  the  Optic  Axis.  Our  third  example  deals  with 
wave  propagation  along  the  optic  axis.  Recall  the  optic  axis  lies  along  the  z-axis.  The 
direction  cosines  describing  wave  propagation  in  the  positive  z-direction  are  given  by 


a  =  0,  a  =  0,  a  =  1. 
x  y  z 


(90) 


Substituting  the  dielectric  tensor  elements  from  equation  (80)  and  the  direction  cosines 
from  equation  (90)  back  into  the  8-by-8  eigenvalue  equation  (64)  gives  the  following 
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eigenvalue  equation: 


0 

0 

0 

0 

0 

-e-1 

0 

0 

0 

zE 

z'E 

X 

0 

0 

0 

0 

e-1 

0 

0 

0 

0 

zE 

y 

z'E 

y 

0 

0 

0 

-e-1 

e 

0 

0 

0 

0 

z'E 

z 

z'E 

z 

0 

0 

e 

e 

0 

0 

0 

0 

0 

0 

2 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

H 

X 

n 

H 

X 

1 

0 

0 

0 

0 

0 

0 

0 

H 

y 

H 

y 

0 

0 

0 

0 

0 

0 

0 

-1 

H 

z 

H 

z 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

(91) 


It  can  be  easily  shown  that  only  one  physically  acceptable  eigenvalue  satisfies  equation 
(91).  The  corresponding  index  of  refraction  is  given  by  the  following: 


The  corresponding  eigenvector  solution  is  given  by 
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where  C  and  C  are  arbitrary  constants.  This  eigenvector  solution  describes 
monochromatic  electromagnetic  plane-waves  propagating  along  the  z-direetion  (optic 
axis)  with  speed  of  cjn.  The  electric  and  magnetic  field  vectors  of  this  wave  are 
orthogonal  to  each  other  as  well  as  the  direction  of  wave  propagation.  Any  state  of 
polarization  is  allowed  providing  it  satisfies  the  transversality  conditions  summarized  in 
equation  (93). 


4.2  OPTICALLY  ACTIVE  MATERIALS 


Certain  optical  materials  are  found  to  possess  the  ability  to  rotate  the  plane  of 
polarization  of  electromagnetic  radiation  passing  through  them.  This  phenomenon  is 
commonly  referred  to  as  optical  activity.  Optical  activity  can  be  explained  on  the  basis  of 
the  assumption  that  the  speed  of  propagation  for  left  circularly  polarized  light  in  the 
material  is  different  from  that  of  right  circularly  polarized  light.  If  a  dielectric  material  is 
placed  in  a  static  magnetic  field  and  a  beam  of  linearly  polarized  light  is  sent  through  the 
dielectric  material  in  the  direction  of  the  applied  magnetic  field,  then  a  rotation  of  the 
plane  of  polarization  of  the  emerging  light  is  found  to  occur.  In  other  words,  the  presence 
of  the  magnetic  field  causes  the  dielectric  material  to  become  optically  active.  This  is 
known  as  the  Faraday  effect. 

It  is  a  simple  matter  to  show  that  if  the  dielectric  tensor  has  conjugate  imaginary  off- 
diagonal  elements  (Fowles,  1968),  namely 
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where  e12is  real,  then  the  tensor  describes  an  optically  active  medium.  For  the  case  of 
the  Faraday  effect,  the  off-diagonal  dielectric  tensor  element  is  proportional  to  the 
applied  static  magnetic  field  strength  assumed  to  lie  along  the  z-axis.  The  optic  axis  of 
the  material  also  lies  along  the  direction  of  the  applied  magnetic  field.  Note  the  dielectric 
tensor  (equation  94)  is  identical  to  that  of  a  uniaxial  crystal  aside  from  the  off-diagonal 
elements. 
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Using  equation  (63),  we  find  that  the  optical  impermeability  tensor  has  the  form 


(95) 


Let  us  now  look  at  another  example  to  illustrate  these  ideas. 

Example  4:  Wave  Propagation  Along  the  Optic  Axis.  For  this  case  we  have 

ax  =  °'  ay  =  °'  az  =  L  (96) 

Substituting  the  dielectric  tensor  (equation  94),  the  optical  impermeability  tensor 
(equation  95),  and  the  direction  cosines  (equation  96)  back  into  the  electromagnetic  wave 
equation  matrix  representation  (equation  61a)  gives  the  3-by-3  eigenvalue  equation 


It  can  be  easily  shown  only  two  physically  acceptable  eigenvalues  satisfy  equation  (97). 
The  corresponding  eigenindices  of  refraction  are  given  by  the  following  formula 
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Associated  with  each  value  of  the  index  of  refraction  n  is  a  corresponding  eigenvector 
solution  (eigenstate  of  polarization).  The  eigenvector  solutions  are 
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where  C  is  an  arbitrary  constant.  These  two  eigenvector  solutions  represent  orthogonal 
states  of  circularly  polarized  electromagnetic  radiation. 


4.3  ELECTROOPTICAL  MATERIALS 

In  this  subsection,  we  consider  the  propagation  of  electromagnetic  radiation  in  a 
crystalline  material  in  the  presence  of  an  applied  static  electric  field.  For  certain  types  of 
crystals,  the  application  of  an  electric  field  results  in  a  change  of  the  optical  birefringent 
properties  of  the  crystal.  This  is  referred  to  as  the  electrooptic  effect.  The  electrooptic 
effect  affords  a  convenient  means  of  controlling  the  phase  or  intensity  of  electromagnetic 
radiation.  The  electrooptic  effect  is  used  in  a  number  of  applications,  including  optical 
modulation,  optical  beam  deflection,  and  spectral  tunable  filters  (Yariv  &  Yeh,  1984). 

There  are  primarily  three  electrooptic  effects  described  to  various  degrees  in  the 
literature.  They  are  (1)  the  Pockels  effect,  (2)  the  Kerr  effect,  and  (3)  the  Stark  effect.  The 
Pockels  effect  is  concerned  with  the  alteration  in  the  refractive  properties  in  a 
piezoelectric  crystal  by  the  application  of  a  strong  electric  field.  The  effect  is  proportional 
to  the  first  power  of  the  electric  field.  The  Kerr  effect  deals  with  the  occurrence  of 
birefringence  in  a  transparent  isotropic  medium  when  it  is  placed  in  an  electric  field.  The 
medium  then  behaves  like  a  uniaxial  crystal  with  its  optic  axis  lying  in  the  direction  of 
the  applied  electric  field.  The  effect  is  proportional  to  the  square  of  the  electric  field 
strength.  The  Stark  effect  concerns  the  displacement  and  splitting  of  the  lines  in  atomic 
spectra,  and  the  appearance  of  new  lines,  owing  to  the  influence  of  a  transverse  electric 
field.  For  pedagogical  purposes,  we  will  only  concern  ourselves  with  the  electrooptic 
Pockels  effect.  A  detailed  discussion  on  the  various  electrooptic  effects  can  be  found  in 
the  book  of  selected  reprints  on  electrooptic  devices  by  Kaminow  (1974). 
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A  popular  method  for  determining  the  eigenindices  of  refraction  and  corresponding 
eigenvectors  associated  with  wave  propagation  through  electrooptical  materials  is 
through  the  use  of  the  index  ellipsoid  (or  optical  indicatrix)  describing  the  electrooptic 
crystal  of  interest  (Yariv  &  Yeh,  1984).  For  the  case  of  the  Pockels  effect,  the  mathematical 
form  of  the  index  ellipsoid  equation  is  given  by 

Tlli:c2  +  T,22^2  +  T133z2  +  2T|23^2  +  2ll31ZX  +  2ri12X^  =  l'  (100) 


where 


111 

i/nx 

rll 

r12  r13 

^22 

l/nl 

r21 

r22  r23 

^33 

y 

2 

r31 

r32  r33 

— 

1  /ni 

+ 

^23 

z 

0 

0 

r41 

r42  r43 

% 

r51 

r52  r53 

<N 

P” 

1 

0 

_r61 

r62  r63 

(101) 


The  quantities  T| . .  are  the  elements  of  the  optical  impermeability  tensor,  r .  are  the  linear 
(or  Pockels)  electrooptic  coefficients,  are  the  components  of  the  applied  static  electric 
field  along  the  three  principle  crystallographic  axes,  and  n .  are  the  principal  indices  of 
refraction  of  the  electrooptic  crystal  in  the  absence  of  the  applied  electric  field. 

An  alternative  method  for  determining  the  eigenindices  of  refraction  and 
corresponding  eigenvectors  is  through  the  use  of  either  the  Maxwell  eigenvalue  matrix 
representation  (equation  58)  or  the  electromagnetic  wave  equation  matrix  representation 
(equation  61).  The  dielectric  tensor  [e]  can  be  determined  from  the  optical 
impermeability  tensor  [q]  with  the  help  of  equation  (57),  namely 

[e]  =  [T|]-1.  (102) 

With  the  help  of  equation  (101),  the  optical  impermeability  tensor  can  be  written  in  a 
form  compatible  with  both  equations  (58)  and  (61),  namely 

[TlJ  =  [Tl0]  +  [Aq],  (103) 


30 


where,  by  definition 
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Let  us  look  at  one  final  example,  now  illustrating  wave  propagation  in  an  electrooptical 
medium. 

Example  5:  Bismuth  Silicon  Oxide  (Bi12SiO20):  Bismuth  silicon  oxide  belongs  to  the 
cubic  class  of  crystals  (Yariv  &  Yeh,  1984),  hence 
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In  addition, 
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All  other  Pockels  (linear)  electrooptic  coefficients  are  equal  to  zero  for  bismuth  silicon 
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oxide,  hence  equation  (105)  simplifies  to 
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Suppose,  for  this  example,  that  the  applied  static  electric  field  points  only  along  the  z- 
axis.  That  is. 


*  0. 


Also,  let  us  consider  wave  propagation  along  the  z-axis.  For  this  case  we  have 


(109) 


a  =  0,  a  =  0,  a  =  1. 
x  y  z 


(110) 


Substituting  the  optical  impermeability  tensor  information  (equations  103,  106,  108  and 
109),  dielectric  tensor  (equation  102),  and  the  direction  cosines  (equation  110)  into  the 
electromagnetic  wave  equation  matrix  representation  (equation  61a)  gives  the  following 
3-by-3  eigenvalue  equation: 
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It  can  be  shown  only  two  physically  acceptable  eigenvalues  satisfy  equation  (111).  The 
corresponding  eigenindices  of  refraction  are  given  by  the  following  formula 


Associated  with  each  value  of  the  index  of  refraction  n  is  a  corresponding  eigenvector 
solution  (eigenstate  of  polarization).  The  eigenvector  solutions  are 
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where  C  is  an  arbitrary  constant.  These  two  eigenvector  solutions  represent  orthogonal 
states  of  linearly  polarized  electromagnetic  radiation. 

5.0  SUMMARY  AND  CONCLUSIONS 


From  the  vector  form  of  the  Maxwell  field  equations  (4),  an  8-by-8  covariant  matrix 
formulation  of  Maxwell's  equatior ^  (8)  was  developed.  Tantamount  to  this  formulation 
is  a  skew-Hermitian  spac>.  c.^e  8-by-8  differential  matrix  operator  [M]  defined  in 
equation  (1).  In  addition,  r  .er  landmark  effects  of  classical  electromagnetic  theory  were 
easily  derived  with  the  use  of  the  8-by-8  matrix  operator  [M]:  (a)  the  electromagnetic 
wave  and  charge  continuity  equations  (11);  (b)  definition  of  the  electromagnetic 
potentials  (employing  the  Lorentz  gauge)  in  terms  of  the  electromagnetic  fields,  and  the 
Lorentz  conditions  (equation  17);  (c)  the  wave  equations  for  the  electromagnetic 
potentials  (equation  23);  and  (d)  Poynting's  theorem  (equation  26)  on  energy 
conservation.  These  matrix-based  derivations  completely  avoid  the  need  for  vector 
calculus  and  attendant  special-purpose  identities  involving  multiple  curls  and 
divergences.  The  necessary  mathematical  baggage  is  strongly  reduced  by  the  approach. 
Note  that  the  extension  of  the  matrix  [M]  into  a  square  form,  described  above,  is  essential 
fo  the  above  electromagnetic  derivations. 

The  matrix  form  of  the  Maxwell  field  equations  (8)  were  then  cast  into  the  Maxwell 
eigenvalue  8-by-8  matrix  representation  (equation  51).  This  eigenvalue  matrix 
representation  was  used  to  solve  for  the  eigenindices  of  refraction  and  corresponding 
polarization  eigenstates  for  a  variety  of  wave-propagation  problems  dealing  with  linear 
homogeneous  anisotropic  optical  media  of  infinite  extent  in  the  presence  of 
monochromatic  plane-wave  electromagnetic  fields.  In  a  similar  manner,  the 
electromagnetic  wave  and  continuity  equations  (11)  were  used  to  formulate  a  wave 
equation  eigenvalue  3-by-3  matrix  representation  (equation  61)  for  nonmagnetic 
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materials.  This  eigenvalue  representation  can  also  be  used  to  find  the  eigenindices  of 
refraction  and  polarization  eigenstates  for  monochromatic  plane-wave  electromagnetic 
fields. 

A  variety  of  examples  were  considered  to  illustrate  the  usefulness  of  these  matrix 
eigenvalue  representations  in  solving  wave-propagation  problems.  In  particular,  the 
eigenindices  of  refraction  and  corresponding  polarization  eigenstates  associated  with 
wave  propagation  in  vacuum,  dielectrics,  uniaxial  crystalline  media,  optically  active 
media  and  electrooptical  media  were  easily  obtained.  Off-the-shelf  computer  software 
packages,  like  MATLAB  for  numerical  computations  and  MATHEMATICA  for  symbolic 
manipulations,  are  well  suited  for  finding  these  eigenindices  and  polarization 
eigenstates. 
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